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A search for a narrow baryonic resonance decaying to K^p or K^p is carried out in 
deep inelastic ep scattering with the HI detector at HERA. Such a resonance could be 
/\^ ■ a strange pentaquark 0+, evidence for which has been reported by several experiments. 

j^ . The K^p and K^p invariant mass distributions presented here do not show any significant 

peak in the mass range from threshold up to 1.7 GeV. Mass dependent upper limits on 
o{ep -^ e0+X) X BR{@^ -^ K^^p) are obtained at the 95% confidence level. 
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1 Introduction 

Recently several fixed-target experiments have published evidence for the production of a strange 
pentaquark^ 0+ [1], a hypothetical baryon [2] with a minimal quark content of uudds, observed 
in the decay channels K^n and K^p. This state has been reported with masses in the range of 
1520 to 1540MeV and with a narrow width, consistent with the experimental resolution in most 
of the observations. Evidence for 0+ production has been also obtained in deep inelastic ep 
scattering (DIS) at HERA by the ZEUS experiment [3]. Many non-observations have also been 
reported [1]. The experimental situation is thus controversial and further data are needed to 
establish the existence of this resonance. 

This paper presents a search for the strange pentaquark 0+ using 74pb~^ of deep inelastic 
ep scattering data taken with the HI detector in the years 1996-2000. A narrow resonance is 
searched for in the K^p or K^p decay channel in the mass range from 1.48 to 1.7 GeV and in 
the kinematic range of negative four momentum transfer squared, Q^, from 5 to 100 GeV^ and 
of inelasticity, y, from 0.1 to 0.6. 



2 Experimental Procedure 
2.1 HI Apparatus 

A detailed description of the HI detector can be found in [4]. The following briefly describes 
only those detector components important for the present analysis. 

The tracks from charged particles used in this analysis are reconstructed in the central 
tracker, whose main components are two cylindrical drift chambers, the inner and outer central 
jet chamber (CJCs). The inner and outer CJC are mounted concentrically around the beam-line, 
covering the range of pseudorapidities^ — 1.9 < r] < 1.9 for tracks coming from the nominal 
event vertex. The CJCs lie within a homogeneous magnetic field of 1.15T which allows the 
transverse momentum, p^, of charged particles to be measured. Two additional drift chambers 
complement the CJCs by precisely measuring the z coordinates of track segments and hence 
assist in the determination of polar angles. Two cylindrical multi-wire proportional chambers 
facilitate triggering on tracks. The transverse momentum resolution of the central tracker is 
o{p^)/Pj — 0.005/?^/ GeV ©0.015. Charge misidentification is negligible for particles orig- 
inating from the primary vertex and having transverse momenta in the range relevant to this 
analysis. 

The specific ionisation energy loss of charged particles, dE/dx, is derived from the mean of 
the inverse square-root of the charge collected by all CJC sense wires with a signal above a cer- 
tain threshold. The average resolution for minimum ionising particles is (7{dE /dx) / {dE /dx) ~ 
8% [5]. 



^In this paper particle names are used to refer to both the particle and its antiparticle, unless explicitly stated 
otherwise. 

^The pseudorapidity is given by 7] = — In tan 0/2, where the polar angle 9 is measured with respect to the z axis 
given by the proton beam direction. 



A lead/scintillating-fibre calorimeter (SpaCal) is located in the direction of the electron 
beam and covers the pseudorapidity range —1.39 < rj < —3.64. It contains electromagnetic 
and hadronic sections and is used here to detect the scattered electron in DIS events and to mea- 
sure its energy. A planar drift chamber, positioned in front of the SpaCal, measures the polar 
angle of the scattered electron track originating from the event vertex. The global properties of 
the hadronic final state are reconstructed combining information from the central tracker, the 
SpaCal and the Liquid Argon calorimeter, which surrounds the central tracker. The DIS events 
studied in this paper are triggered on the basis of an energy deposition detected in the SpaCal, 
complemented by signals in the CJCs and the multi-wire proportional chambers. 

The luminosity measurement is based on the Bethe-Heitler process ep -^ epf, where the 
photon is detected in a calorimeter located downstream of the interaction point. 



2.2 Monte Carlo Simulation of 0"^ 

To estimate the acceptance for the detection of a hypothetical 0+ state, a Monte Carlo simula- 
tion based on the RAPGAP 3.1 [6] event generator is used, incorporating the Lund string model 
fragmentation [7] as implemented in PYTHIA 6.2 [8]. The kinematic distributions of strange 
baryons in DIS data are reasonably well described [9] by RAPGAP. The 0+ is introduced by 
changing the mass of the E*+ to values in the range from 1.48 to 1.7 GeV and forcing it to decay 
to K^p. By doing so, it is assumed that the 0+ is produced at pseudorapidities and transverse 
momenta similar to those of other strange baryons and that it decays isotropically. In this sim- 
ulation the 0+ particle is produced on mass shell. The generated events are passed through the 
HI detector simulation based on GEANT [10] and are then subjected to the same reconstruction 
and analysis chain as are the data. 



2.3 Selection of DIS Events 

The analysis is carried out using data corresponding to an integrated luminosity of ^ = 74pb^ ^ , 
taken in the years 1996-2000. During this time HERA collided electrons-^ at an energy of 
27.6 GeV with protons at 820 GeV (1996-1997) and 920 GeV (1998-2000)^. 

Events are selected if the z coordinate of the event vertex, reconstructed using the central 
tracker, lies within 35cm of the mean position for ep interactions. The scattered electron is 
required to be reconstructed in the SpaCal with an energy. Eg, above 1 1 GeV. The negative four 
momentum transfer squared of the exchanged virtual photon, Q^, is required to lie in the range 
5 < Q^ < 100 GeV^, as reconstructed from the energy and polar angle of the scattered electron. 
The inelasticity y of the event is reconstructed using the scattered electron kinematics and is 
required to be in the range 0.1 < y < 0.6. The lower cut on y ensures that the hadronic final 
state lies in the central region of the detector, whilst the upper cut corresponds approximately 
to the cut on Eg. The difference between the total energy E and the longitudinal component of 

^The analysis uses data from periods when the beam lepton was either a positron (^ — 65 pb^') or an electron 
(^ = 9pb-i). 

^The sample with a proton energy of 820 (920) GeV corresponds to a luminosity of ^ =18 (56) pb^^ resulting 
in an effective v^ = 314 GeV for the total sample. 



the total momentum p^, calculated from the electron and the hadronic final state, is restricted 
io35 <E — p^<10 GeV. This requirement suppresses photoproduction background, in which 
the electron escapes detection and a hadron fakes the electron signature. 



'0 



2.4 Selection of K, Meson and Proton Candidates 
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The analysis is based on charged particles reconstructed in the central tracker. Tracks are ac- 
cepted if they have transverse momenta/?^ > 0.15 GeV and pseudorapidities \r\\ < 1.75. The 
K^ meson is identified through its decay into charged pions, K^ -^ tz^tz^ . Events are accepted 
if they contain at least one K^ candidate and at least one proton candidate track originating from 
the primary vertex. 

K^ candidates are searched for by performing a constrained fit to each pair of oppositely 
charged tracks. The fit demands these tracks to originate at a common secondary decay vertex 
and the decaying neutral particle to come from the primary vertex. The secondary vertex must 
be radially displaced by at least 2 cm from the primary interaction point. The candidates are re- 
quired to have a transverse momentum p^ {K^) > 0.3 GeV and a pseudorapidity 1 7] {K^) \ < 1.5. 
Contamination from A production is eliminated by requiring that the invariant mass Mp^^ of 
the two tracks, reconstructed assigning the proton (pion) mass to the track with higher (lower) 
momentum, be above 1.125 GeV. Background from converted photons is rejected by the re- 
quirement Mee > 50MeV. Figure [It shows the distribution of the invariant mass Mj^+j^- of 
the K^ candidates together with a fit to the data using a superposition of two Gaussian func- 
tions (to account for different invariant mass resolutions in different decay topologies) and a 
straight line to approximate the background. The fitted peak position is M^+^- = 495.9 MeV 
which agrees with the nominal Kf mass [11] within a few per mill. 133,000 K^ candidates 
are reconstructed, as given by subtracting the fitted background from the data. K^ candidates 
with 0.475 < Mj^+j^- < 0.515 GeV are selected for further analysis. In this mass range the 
background under the K^ peak is ~ 3%. 

Proton candidates are selected using requirements on the specific ionisation energy loss, 
dE/dx , measured in the CJCs. Figure [1J5 shows the measured dE/dx plotted against momen- 
tum for all tracks originating from the primary vertex, which lead to a mass M^o„ < 1.8 GeV 
when combined with the K^ candidates. The curves in Fig. [ij? represent the most probable 
dE/dx values as derived from a phenomenological parameterisation [5] based on the Bethe- 
Bloch formula. Likelihoods for a particle to be a pion, kaon or proton are obtained from the 
difference between the measured dE/dx and the most probable value for each particle type at 
the reconstructed momentum. The normalised proton likelihood, Lp, is defined as the ratio of 
the proton likelihood to the sum of the pion, kaon and proton likelihoods. In order to optimise 
simultaneously the background suppression and the proton selection efficiency, a momentum 
dependent cut on the normalised proton likelihood Lp is applied of L^ > 0.3 {Lp > 0.1) for 
proton momenta below (above) 2 GeV. The efficiency of the dE/dx selection is tested using 
protons from A decays. The efficiency varies between 65% and 100% as a function of momen- 
tum and is described by the Monte Carlo simulation to within 5%. 
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Figure 1 : a) Inclusive K^ signal in the invariant n^n^ mass distribution for 5 < Q^ < 100 GeV^ 
together with the result from a fit of a sum of two Gaussian functions for the signal and a straight 
line for the background, b) Specific ionisation energy loss relative to that of a minimally ionising 
particle plotted as a function of momentum. The lines indicate parameterisations of the most 
probable energy loss for pions, kaons and protons measured in the CJCs. 

3 Analysis of K^p Combinations 

In order to search for a 0+ resonance, the candidate K^ mesons are combined with the pro- 
ton candidates. To improve the mass resolution the K^p four vector is calculated as the sum 

of the K^ and proton four vectors with E^o = ^ /p^o +^^0' where the nominal mass M^o is 

used instead of Mj^+j^-. For the K^p system, p^{K^p) > 0.5 GeV and \r\{K^p)\ < 1.5 are re- 
quired. The M^o distributions are shown in Fig. |2lfor three bins in Q^ (5 < Q^ < 10 GeV^, 
10 < 2^ < 20 GeV^ and 20 < Q^ < 100 GeV^). The shape of the invariant mass distributions 
is found to be reproduced by a background Monte Carlo simulation of inclusive DIS events 
using the DJANGOH event generator [12] and the HI detector simulation based on GEANT. A 
fit of an empirical background function of the form 



/3 



fiMKOp) = a • (Mk^p -MthrY ■ exp{-(M^Op -M,/,,) ■ 7} 



(1) 



is performed to the data, where Mthr = Mj^o +Mp (Mp being the proton mass) and a, j8 and 7 

are free parameters determined for each Q^ interval independently. The data are well described 
by this phenomenological function. No narrow resonance is observed in any of the Q^ bins. 
The Mf^ distribution is therefore used to set upper limits on the 0+ production cross section. 



C7^ = {a{ep -^ e0+X) +a{ep -^ e@+X)) x 5i?(0+ -^ K^p) . 

Since the mass of the 0+ candidate is experimentally not well established, mass dependent 
limits are derived in the range from 1.48 to 1.7 GeV. For a given 0+ mass, M@+, the expected 
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Figure 2: Invariant K^p{p) mass spectra in bins of Q^. The full line shows the result from 
the fit of the background function ([T]) to the data. The upper limits on the cross section C7®^ 
(see text) at 95% confidence level integrated over the kinematic range p^{K^p) > 0.5 GeV, 
1 7] (K^p) I < 1 .5 and 0.1 < y < 0.6 are shown below the mass spectra. 



number of selected K^p combinations due to 0+ production is related to o® via 

N{Mq+)=o®-^-£dis-£&+{Mq+)-BR{K^^k+k-) , (2) 

where ^ is the integrated luminosity, Edis is the acceptance of the inclusive DIS event selection 
and £@+{Mq+) is the acceptance of the 0+ selection. The cross section a® is integrated over 
the visible kinematic range studied, which is given by p^{K^p) > 0.5 GeV, \ri{K^p)\ < 1.5, 
0. 1 <y <0.6 and the respective Q^ bin. The branching ratio for the transition of K^ to K^ and 
its decay into charged pions is BR{K^ -^ n+n-) = BR{K^ ^ K^) x BR{K^ ^ n+n-) = 
0.5x(0.6895±0.0014)[ll]. 

An upper limit at the 95% confidence level (C.L.) on N{Mq+), Nul{M@+), is obtained from 
the observed, the background and the signal M^o„ distributions in the mass range from 1 .45 
to 1.8 GeV, using a modified frequentist approach based on likelihood ratios [13]. This takes 
into account statistical and systematic uncertainties of the signal and the background number 
of ^^^;? combinations. The Mj^o„ distribution for signal combinations is taken to be a Gaussian 
with a mean Mq+ and a width corresponding to the experimental mass resolution as obtained in 
the 0+ Monte Carlo simulations. This width o{Mq+) varies from 4.8 to 1 1.3MeV in the mass 
range from 1.48 to 1.7 GeV. The background M^o„ distribution is taken to be the fitted function 
given by equation (U]). 

A systematic uncertainty on this background distribution is assessed by performing the fit 
under different assumptions: using the background function ([T]) in the full mass range, excluding 
a mass window of ±2c7 around the 0+ mass, and also using the sum of the background function 
and a Gaussian with fixed mass Mq+ and width o{M@+) to account for a possible signal. The 
uncertainty of the number of background K^p combinations is estimated from the difference 
between the different fitting methods and amounts to 2%. 

The systematic uncertainty of N{Mq+) comprises the following main contributions: 

• The measurement of the luminosity has an uncertainty of 1.5%. 

• The uncertainty of the inclusive DIS event selection, Edis^ is 6.5%, which is coming 
mainly from contributions due to the trigger efficiency (5%), the SpaCal energy cali- 
bration (3%), remaining contamination from photoproduction background (2.5%) and 
radiative corrections (1%). 

• The efficiency of the 0+ selection, £0+ {Mq+ ) , has an uncertainty of 8% which comprises 
the uncertainty in modelling track losses (6%) and the uncertainty in the efficiency of the 
dE/dx selection (5%). 

• The Monte Carlo model used for correction is based on the assumption that pentaquarks 
are produced with similar phase space distributions as strange baryons. Since no es- 
tablished production mechanism for the 0+ yet is known, production model dependent 
uncertainties are not considered. Dependencies on the QCD models are estimated by 
comparing the 0+ acceptances derived with the RAPGAP and the CASCADE [14] event 
generators, which incorporate different QCD evolution schemes. The difference is found 
to be small and negligible compared with other sources of systematic uncertainties. 



The contributions are added in quadrature and the resulting total systematic uncertainty 
of N{Mq+) is 11%. The uncertainty of the number of background K^p combinations is the 
dominating systematic effect in the limit calculation. 

The upper limit on the cross section, a®^, is then calculated from N[jl{Mq+) according to 
equation ©. The upper limits at 95% C.L. are shown below the mass spectra of Fig.|3for the 
three different Q^ bins. The limits vary between 30 and 90 pb for the different Q^ bins and 
over the mass range from 1 .48 to 1 .7 GeV. The invariant mass spectra of positive K^p and 
negative K^p combinations are also studied separately. No narrow resonance is observed. The 
corresponding upper limits for the 0+ decaying to K^p and its charge conjugate 0+ decaying 
to K^p, shown in Fig.|3l are found to be of comparable size. The up- and downward fluctuations 
of the limits occur at different masses for the different Q^ bins and charges, which supports the 
hypothesis that the observed K^p invariant mass distributions are consistent with being due to 
combinatorial background only. 

The ZEUS experiment has reported a positive 0+ observation at a mass of 1.522 GeV in 
DIS for Q^ > 20 GeV^ using a data sample corresponding to an integrated luminosity of 121 
pb^^ [3]. The K^p system was reconstructed using only protons having a momentum below 
1.5 GeV, while the requirements on the transverse momenta and pseudorapidities of ^^ and 
K^p are the same as in the present analysis. The analysis described above is repeated using only 
proton candidates with momenta below 1.5 GeV. The resulting invariant K^p{p) mass spectra 
are shown in Fig. |4^ for 20 < 2^ < 100 GeV^ and 0.1 < y < 0.6. No significant pentaquark 
signal is observed for events in the low momentum proton selection. The upper limits a®^ at 
95% C.L., derived from these mass spectra, are shown in Fig.|4j5. At a 0+ mass of 1.52 GeV an 
upper limit on the cross section of 72pb at 95% C.L. is found. The K^p and K^p combinations 
do not yield any significant peak either. The corresponding upper limits are also shown in 
Fig. St. 



4 Conclusions 

A search for the strange pentaquark 0+ in deep inelastic ep scattering is presented. No signal 
for 0+ production is observed in the decay mode 0+ -^ K^p and 0+ -^ K^p for negative 
momentum transfers squared, Q^, between 5 and 100 GeV^. Assuming that pentaquarks are 
produced with similar kinematics as known strange baryons, mass dependent upper limits at 
95% confidence level on the cross section a{ep -^ e&^X) x BR{&^ -^ K^p) are derived in 
intervals of Q^ and found to vary between 30 and 90 pb over the mass range from 1.48 to 
1.7 GeV. 

The analysis is repeated, restricted to large Q^ and low proton momentum, a region in which 
the ZEUS collaboration observes evidence for a 0+ signal. For this selection no signal is found 
either. 
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Figure 3: Upper limits on the cross section c®^ (see text) at 95% confidence level in bins of 
Q^ for K^p (full line) and K^p (dashed line) separately, integrated over the kinematic range 
Pj{K^p) > 0.5 GeV, \r]iK^p)\ < 1.5 andO.l < j < 0.6. 
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Figure 4: a) Invariant K^p{p) mass spectra for 20 < Q^ < 100 GeV^ for proton candidates with 
momenta below 1.5 GeV and b) upper limits on the cross section C7®^ (see text) at 95% confi- 
dence level for 20 < Q^ < 100 GeV^ integrated over the kinematic range p^{K^p) > 0.5 GeV, 
\r\{K^p)\ < 1.5 and 0.1 <y <0.6 using the low momentum proton selection for all K^p com- 
binations and c) for K^p (full line) and K^p (dashed line) combinations separately. 
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